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We consider a one-dimensional quantum many-body system and investigate how the interplay 
between interaction and on-site disorder affects spatial localization and quantum correlations. The 
hopping amplitude is kept constant. To measure localization, we use the number of principal com- 
ponents (NPC), which quantifies the spreading of the system eigenstates over vectors of a given basis 
set. Quantum correlations are determined by a global entanglement measure Q, which quantifies 
the degree of entanglement of multipartite pure states. Our studies apply analogously to a one- 
dimensional system of interacting spinless fermions, hard-core bosons, or yet to an XXZ Heisenberg 
spin-1/2 chain. Disorder is characterized by both: uncorrelated and long-range correlated random 
on-site energies. Dilute and half-filled chains are analyzed. In half-filled clean chains, delocaliza- 
f"*^ . tion is maximum when the particles do not interact, whereas multi-partite entanglement is largest 

when they do. In the presence of uncorrelated disorder, NPC and Q show a non-trivial behavior 
with interaction, peaking in the chaotic region. The inclusion of correlated disorder may further 
extend two-particle states, but the effect decreases with the number of particles and strength of 
their interactions. In half-filled chains with large interaction, correlated disorder may even enhance 
C"j | localization. 
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I. INTRODUCTION 



Disorder may significantly affect the properties of physical systems. Spatial localization of one-particle states 
(Anderson localization), for example, is due to uncorrelated random disorder [J 0, H, 0| ; whereas short range @, [E0 
and long range correlated @, || disorder have been associated with the appearance of delocalized states. Interestingly, 
however, a recent experiment [10 ] shows that long-range correlated disorder may in fact either suppress or enhance 
localization. This scenario becomes yet more complex when two or more particles are considered and the effects of 
interactions are taken into account. A clear picture of the interplay between interaction and disorder is essential 
to advance our understanding of thermodynamic and kinetic properties, as well as the dynamical behavior of real 
\^ quantum many-body systems. 

I/* - ) . The interest in how interactions influence localization and transport properties of mesoscopic systems was boosted 
t^J- ■ with the realization that persistent currents could only be explained if the role of the electron-electron interaction 
was addressed fllTl . In general, interaction between particles in disordered systems may lead to delocalization [l2l. [T3. 
HH EE EE EE 13 EE EE HII ■ The competition between interaction and disorder in many-body systems relates also 
to other interestingphenomena, such as the rich variety of quantum phase transitions of ultracold atomic Bose gases 
in optical lattices I22 |, the transition from integrability to chaos [H,[24 J25[ in spin systems, and the enhancement of 
entanglement [26l 1271 l28l |29| ] and the 'melting' of quantum computers 30j in the context of quantum information. 
• rH , In the present work we investigate how disorder combined with interactions affect spatial localization and quantum 
' correlations of a one-dimensional quantum system. In addition to uncorrelated disorder, we analyze also long-range 
correlated disorder. Correlated disorder may appear in real systems [31] and may also be engineered. The latter 
includes the introduction of scatterers fiol or s pec kles [32| in a one-dimensional waveguide or yet the individual 



tuning of on-site energies via local fields [33j, |34|, |35|. In contrast to the extensively studied Hubbard model, where 
interaction occurs between particles in the same site, we consider interaction between particles in neighboring sites. 
Both dilute and half-filled chains are analyzed. We study the level of delocalization, determined by the number 
of principal components NPC, and the amount of multi-partite entanglement, quantified by a global entanglement 
measure Q, of all eigenvectors of the system. Our numerical results show that: (i) In a non-interacting system, NPC 
and Q decrease with uncorrelated on-site disorder and increase with correlated on-site disorder; (ii) In a clean system, 
interaction restricts delocalized states to narrow energy bands and the average spatial delocalization decays, while 
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the behavior of quantum correlations is non-monotonic; (iii) In a disordered system the behavior of NPC and Q with 
interaction depends highly on the number of particles, being non-trivial in a half- filled chain, where quantum chaos 
may develop; (iv) On average, states of two-interacting particles appear to be delocalized in the considered finite 
systems with uncorrelated disorder, while one-particle states localize. 

The paper is organized as follows. Sec. II describes the model, the relation determining on-site disorder, and the 
quantities computed. The half-filled chain is studied in Sec. Ill, disordered systems with non-interacting, weakly 
interacting and strongly interacting particles are compared. Sec. IV considers the dilute limit and compares the 
results for one and two-particle states. Discussions and concluding remarks are presented in Sec. V. 



II. SYSTEM MODEL AND MEASURES OF DELOCALIZATION AND QUANTUM CORRELATIONS 



A. System Model 



The analysis developed here applies to different one-dimensional quantum many-body systems, including XXZ 
Heisenberg spin-1/2 models and chains of interacting spinlcss fermions or hard-core bosons. 

The Heisenberg spin-1/2 chain describes quasi-one-dimensional magnetic compounds [36| and Josephson-junction- 
arrays [37t |38|. It has also been broadl y us ed as a model for proposals of quan tum computers, including those based 
on semiconductor quantum dot arrays [39| , solid state NMR spin systems |40l ] , and electrons floating on Helium [H, 
[34j . We investigate a chain with open boundary conditions and nearest-neighbor interactions, as determined by the 
Hamiltonian: 



H — H i + Hxy — Hq + H int + Hxy, 

L 



71 = 1 

L-l 



Hint JAS^Sj l+1 , 
n=l 

L-l 

Hxy — 2^ J (^n^n+l + ^n^n+l) 



(1) 



Above, h is set equal to 1, L is the number of sites, and S n = a n /2 is the spin operator at site n, a*' y,z being the Pauli 
operators. The parameter Q n = uj + uj n , where lu h = de n , corresponds to the Zeeman splitting of spin n, as determined 
by a static magnetic field in the z direction. In a clean system, all sites have the same energy splitting (d = 0) , whereas 
disorder is characterized by the presence of on-site defects (d ^ 0). The relation specifying e„ is discussed in Sec. II. B 
- correlated and uncorrelated random disorder are considered. J is the exchange coupling strength and A is the 
anisotropy associated with the Ising interaction S*S*+i- We set J, A > 0. The total spin operator in the z direction, 
S z = J2n=i &n> ^ s conserved, therefore the matrix H is composed of independent blocks. Each block belongs to a 
single S z subspace characterized by a fixed number M of total spins pointing up, the dimension N of each sector is 
given by N = (£) = Ll/[(L - M)\M\ _ 

The model of Eq. [T] may be mapped into a spinless fermion system via a Jordan- Wigner transformation [4l| , so 
that Hxy becomes #hop and 

H = H 0i + Hhop = Hq + Hint + -ffhop, 
L 

H = ^ ^n a ™ a " 

71=1 

L-l 

Hint = X! J ^ a li a li+l a n+l a 7i, 

71=1 

i- 1 J 

H h op = 2" (»n a n+l + a l+l a «) > ( 2 ) 

71 = 1 
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Above, ajj and a n are creation and annihilation operators, respectively. The presence of a fermion on site n corresponds 
to an excited spin, the on-site fermion energies are the Zeeman energies, J is the fermion hopping integral, and J A 
gives the fermion interaction strength. The conservation of S z translates here into conservation of the total number 
of particles. Disordered systems of interacting fermions simplified by neglecting the spin degrees of freedom has been 
vastly considered as a first approximation in studies of the metal- insulator transition [15l . |42j • 

The fermionic model above may also be mapped onto a system of hard-core bosons. The Hamiltonian has the same 
form of Eq. [5J where the fermionic creation and annihilation operators are substituted by bosonic ones [43| . A system of 
hard-core bosons in one-dimensional optical lattices constitutes a versatile tool in studies of various complex quantum- 
physical phenomena, and has been receiving enormous attention, specially after experimental realizations [441 ]. 

Here, a particle or an excitation will generically refer to an excited spin, a spinless fermion, or a hard-core boson, 
depending on the specific system one addresses. We assume L even and study both a half- filled (M — L/2) and a 
dilute (M = 2) chain. 



B. On-site disorder 



There are widely different physical, biological, and economical processes modeled as stochastic time series. These 
series may be generated, for example, by imposing a power law power spectrum S(f) oc f~ a 45, 46]. The value of 
a determines the type of noise: a = 0, 1 and 2 correspond, respectively, to white noise (flat frequency spectrum), 
1/f noise, and Brownian noise. We borrow from these studies the specific sequence of long-range correlated on-site 
energies considered in this work, 



L/2 

e» = 

k=l 



2tt 

T 



/ 2nnk 
cos I — - h (pk 



(3) 



where 4>k are random numbers uniformly distributed in the range [0,27r]. The sequence is constructed so that, by 
Fourier transforming the two-point correlation function (e„e m ) , one obtains the power law spectral density S(k) ex k~ a . 
When a — 0, e„'s in w n — de n are random numbers with a Gaussian distribution, leading to the scenario of uncorrelated 
disorder: (uj n ) = and (uo n Ld m ) = d 2 5 n . m . Correlated on-site energies appear for a > 0. The energy sequence is 



normalized: (e„) = and the unbiased dispersion y X)n=i( e « — ( e n)) 2 /(L — 1) = 1. 

Long-range correlations are widespread in biological physics and have been extensively analyzed in this context 
[47l ]. In the field of condensed-matter physics, most works dealing with short-range d, H, 0| and long-range i_H Q 
correlated disorder in the context of mobility edges have been limited to the case of a single particle. Ref. for 
instance, considered sequence ([3]) and showed that when a > 2 the one-particle wave functions remain delocalized 
even in the thermodynamic limit. Here, two or more excitations are taken into account and we investigate how the 
disorder parameters d and a and the interaction amplitude A affect derealization and multipartite entanglement in 
the finite systems described above. 



C. Derealization 



To quantify the extent of derealization of an eigenvector — J2k=i cl j\f k ) °f Hamiltonian (Ill2p written in the 
basis \ip k ), we consider the number of principal components (NPC) [48j|, defined as 



NPC 7 = . (4) 

This quantity is also commonly referred to as inverse participation ratio. 

A large NPCj is associated with a delocalized state where many basis vectors give a significant contribution to the 
superposition whereas a small NPCj is related to a localized state. Clearly, the components of the eigenvectors 
depend entirely on the choice of basis in which to express them. Our approach here is to take a physically motivated 
basis in order to study Anderson localization in disordered systems with interacting particles. Anderson localization 
refers to the exponential spatial localization of wavefunctions. This justifies our choice to consider the site basis, 
which corresponds to a basis consisting of the eigenstates of Hq. In this basis, the interaction H lnt contributes to the 
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diagonal elements of the Hamiltonian, whereas H nop constitutes the off-diagonal elements, the latter being responsible 
for transfer of excitations along the chain. 

Because of its basis-dependence, NPC is not an intrinsic indicator of quantum chaos. To characterize the onset 
of chaos, we consider the level spacing distribution (see Sec. II. D). In studies of atomic and nuclear phy sics, the 
basis-dependence of quantities to measure the complexity of wavefunctions has long been pointed out [49(. It has 
been argued that the mean-field basis is the preferred representation [lij [H(J [EH , separating global properties from 
local fluctuations and correlations of the wavefunctions. Therefore in Sec. III. A. 2, we also briefly compare the results 
for NPC in the site basis with those for two other representations: free particles (FP)-basis and interacting particles 
(IP)-basis. The first consists of the eigenstates of -f/hop _ the integrable Hamiltonian describing a clean system of 
free particles - and the second consists of the eigenstates of -f/hop + -Hint ~ the Hamiltonian for a clean chain with 
interacting particles. -ffh op + Hint is a l so an integrable model and is solved with the Bethe Ansatz method (52j. 

Maximum derealization, NPC ~ N/3, where N is the matrix dimension [ll|, [H, HH , is obtained with states from 
a chaotic system described by random matrices of a Gaussian Orthogonal Ensemble (GOE). The hamiltonian studied 
here may also lead to chaos, but it is a banded random matrix, having only two-body interactions. In this case, the 
large values of NPC may be reached only in the middle of the spectrum, the borders showing less delocalized states, 
as typical of Two-Body Random Ensembles (TBRE) [H, [Hj]. In TBRE's, the local density of states (LDOS) as a 
function of energy is Gaussian and peaked at the center of the spectrum. 

D. Quantum Chaos 

In quantum systems, integrable and non-integrable regimes may be identified by analyzing the distribution of 
spacings s between neighboring energy levels [Til . [5r| . Quantum levels of integrable systems tend to cluster and 
are not prohibited from crossing, in this case the typical distribution is Poissonian Pp{s) = exp(— s). In contrast, 
chaotic systems show levels that are correlated and crossings are strongly resisted, here the level statistics is given 
by the Wigner-Dyson distribution. The exact form of the distribution depends on the symmetry properties of the 
Hamiltonian. In the case of systems with time reversal invariance [57j |. it is given by Pwd(s) = its/2 exp(— its 2 /A). 

For the purpose of illustration, instead of showing the level spacing distributions for the whole range of parameters 
analyzed here, we compute P(s) and associate with it a number given by the quantity rj, which is defined as 

_ J°°lP(s)-P WD ( S )}d S 

71 ~ f;°[p P (s)-p W D(s)}ds> ( > 

where sq ~ 0.4729 is the first intersection point of Pp and Pwd- This quantity, introduced in [58j |. simplifies the 
visualization of the transition from integrability to chaos. For an integrable system: rj — > 1, while for a chaotic system: 
i] — ► 0. In Ref. [59|, rj = 0.3 was taken as an arbitrary value below which the s yste m may be considered chaotic. To 
derive meaningful level spacing distributions, besides unfolding the spectrum [TTl [5r| . all trivial symmetries of the 
system need to be identified. The distributions are computed separately in each symmetry sector. 

In addition to the conservation of S z , the model described by Eq. |T]) in the absence of disorder may also exhibit 
the following symmetries [28| : invariance under lattice reflection, which leads to parity conservation, and, when the 
system is isotropic, conservation of total spin S 2 = (X)«=i Sn) 2 , that is, [H, S 2 ] — (S 2 symmetry). Notice that only 
reflection symmetry exists in a chain with open boundary conditions, whereas reflection and translational symmetries 
would occur in a ring. 

E. Quantum correlations 

To quantify global quantum correlations, we consider the so-called global entanglement, as proposed by Meyer and 
Wallach [6(|. This is a multi-partite entanglement measure employed for lattices of two level systems (qubits). For a 
pure state \ipj) of a chain with L qubits, it is defined as 

2 L 

Q J =2--Y J Tr{p 2 n ), (6) 

n=l 

where p n stands for the density matrix of the chain after tracing over all qubits but n. Qj is therefore linearly related 
to the average purity of each two-level system, that is, it is an average over the entanglements of each qubit with the 
rest of the system [6y, HH, [Hj] . Maximum global entanglement corresponds to Q = 1 . 
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Given the conservation of S z in H ([1} [of the number of particles in H <j2j) ] , this expression may be further simplified. 
For spins: 

Q, = i-^El(^-Kr^}| 2 - (7) 

n=l 

For spinless fermions: 

1 L 

Q^ = l--^|(^|2ata n -l|^)| 2 . (8) 

n=l 

A more general entanglement measure, the so called generalized entanglement (GE), was proposed in [6^ |. GE is 
based on the relationship of the state with a distinguished set of observables, rather than a distinguished subsystem 
decomposition. It coincides with 1 — Qj when the observable set corresponds to the expectation values of the local 
magnetizations, \(ibAo' s!\il>4 ) with /3 = x,y, z. Global entanglement is also closely related with the more broadly used 
von Neumann entropy [281 ]. All of the above are measures of multipartite entanglement, which differ significantly from 
measures, such as concurrence [64| . which aim at capturing pairwise correlations. 

The meaning of global entanglement is probably better understood with examples. Consider, for instance, a 
bipartite system consisting of two qubits. Each qubit may be in state |0) or |1). States such as |00), |01), or even 
[| 00) + 1 01} + |10) + |ll)]/2 are not entangled, Q = 0. The latter, in particular, is maximally delocalized, but not 
entangled, since it may be written as a product state [|0) + |l)]/\/2® [|0) + |l)]/\/2. In contrast, a state such as 
ai|01) +a2|10), where |ai| 2 , \a2\ 2 ^ 0, cannot be written as a product of states of the qubits; the qubits are then non- 
locally correlated and the state is entangled, Q ^ 0. Maximum entanglement, Q = 1, occurs when |ai| 2 = |a.2 1 2 = 1/2, 
which corresponds to the so-called EPR or Bell state. In this case, Tr{p\) — 1/2, that is, the reduced state of each 
qubit is maximally mixed. But global entanglement goes beyond bipartite entanglement and quantifies the amount of 
multi-partite entanglement. The GHZ state, [|00 ... 00) + 111 . . . H)]/\/2, or state [|00> + |11)]/V5® [j00> + |H)]/a/2 
are examples of states with maximum global correlation [611 ] . 

III. HALF-FILLED CHAIN 

We consider a half-filled one-dimensional system with L = 12 and M = 6. This choice corresponds to the largest 
subspace of the Hamiltonian, the sector where chaos sets in first. 

A. Uncorrelated random disorder 

1. Site-basis 

We start by investigating how uncorrelated disorder affects the spatial localization of systems with interacting 
particles. In the main panel of Fig. [T] we plot the average NPC in the site-basis vs. the amplitude d/J of uncorrelated 
Gaussian disorder for different values of the interaction strength. The maximum spatial derealization occurs in a 
clean system in the absence of interaction. At d/J = 0, the level of derealization decreases with A, and for A = 0, 
NPC also decreases with d/J. On the other hand, when both interaction and disorder are present, the behavior is 
not monotonic and NPC reaches a peak for d/J < 1, eventually decreasing again as the disorder becomes larger than 
the hopping integral. 

At d/J = 0, the wavefunctions of the chain with interacting particles are delocalized in the site-basis, although 
the system is integrable. The addition of disorder in cases where < J A < J further increases NPC by breaking 
symmetries and allowing for couplings between more basis states. This is completely antagonic to the behavior 
of chains with non- interacting particles or with very strong interactions (JA 3> J), where disorder only localizes 
wavefunctions. By comparing the curves of NPC with rj in panels A and B, one sees that the derealization peak 
is directly related to the onset of quantum chaos. There, the minimum value of rj for A = 0.5, 1 and 1.5 occurs at 
d ~ J/4. However, it is important to emphasize that the bump in NPC should not be taken alone as an indication 
of chaos. In fact, a similar but less dramatic behavior is observed for the NPC of a two dimensional system, which is 
chaotic already at d/J = (28|. 

To better understand the further derealization of wavefunctions in interacting ID systems when disorder is con- 
sidered, we select the case where A = 2 and compare, for three values of d/J, the histograms of the diagonal elements 
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FIG. 1: Interplay between interaction and uncorrelated random disorder. Main panel: NPC is the site-basis vs. d/ J . Average 
over 924 states and 20 realizations. Panels (A) and (B): Level of chaoticity vs. d/J (unfolding as described in 25, 65]); average 
over 20 realizations. Panels (a;,,-), (b;, r ), (Q, r ): A = 2; left: normalized histogram for the diagonal matrix elements Hoi (Nt is 
the number of states); right: NPC vs. eigenvalues; 20 realizations. All panels: L = 12, M = L/2, Gaussian random numbers: 
a = 0. 



-Hoi in panels a;, b;, c; and the plots of the level of derealization of the eigenvectors of Hamiltonian (jT} [2J vs. their 
corresponding energies E in panels a r , by, and c r . 

When d/J = 0, the energies of the basis vectors form narrow bands of resonant states (aj), the energies being 
determined by the number of pairs of neighboring excitations and by the number of excitations placed at the edges 
of the chain (border effects). These energies range from — JA(L — l)/4 to JA(L — 3)/4 and the bands are separated 
by J A. As the hopping term is turned on, only basis vectors belonging to the same symmetry sector can couple to 
form the wavefunctions. In addition, if J A > J, the effects of J on states belonging to different bands are negligible 
and the bands, although now of finite width, remain separated. In this case, only states placed in the same band 
and showing the same symmetries can mix, so the level of derealization is reduced to the number of states in the 
band. Moreover, the number of states in a band decreases with the size of the clusters, that is, the number of 
neighboring excitations fljj . The overall effect in clean systems is therefore the monotonic decrease of average NPC 
with interaction. The right panel (a r ) shows the values of NPC. For this integrable system, no clear relationship 
between NPC and E is seen, as expected due to the absence of level repulsion. It is interesting to contrast this 
behavior with the 2D clean system, which is chaotic and outlines the shape typical of chaotic TBRE [j5j|] already at 
e?/J = (data not shown). 

By slightly increasing d/J, the bands broaden. If J A > J, they remain uncoupled (b;): the number of resonant 
intra-band states then decreases and so does NPC (b r ). This explains the valley that precedes the peak in the NPC 
curves shown in the main panel. Larger disorder is then needed to overlap the bands (c;) and increase derealization 
(c r ). Notice that the average NPC for d/ J = and 0.3 is approximately the same, although the NPC dependence on 
energy is significantly different (cf. a r and c r ). The integrable clean chain shows no relationship between NPC and 
E (a r ), but as complexity increases and 77 decreases, a relationship similar to those appearing in TBRE's becomes 
evident (c r ): delocalized states are in the middle of the spectrum and only localized states appear in the edges. 

The case of A = 2 is at the borderline, where band overlapping is significant and where level spacing distributions 
close to a Wigner-Dyson distribution may still be obtained. For larger interactions, such as A = 5, the bands are 
very separated at d/J = and do not merge together by increasing d/J, instead, larger disorder simply prevents 
resonances leading to the monotonic decay of NPC (see main panel) . 

Fig. [5] shows the behavior of global entanglement versus uncorrelated disorder for various values of the anisotropy 
[66|. A clean isotropic chain has Q = 1, since the Hamiltonian is invariant under a global rotation of 180 degrees 
around the x-axis [2g|. For <f/J — > 0, contrary to NPC, the largest value of Q is found for a clean system in the 
presence of weak interaction, < A < 1. This suggests that interactions are key ingredients for the generation of 
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FIG. 2: Interplay between interaction and uncorrelated random disorder. Average global entanglement vs. d/J. Average over 
924 states and 20 realizations; L = 12, M = L/2; a = 0. 

non-local correlations in a half-filled chain. Therefore, we identify two main factors contributing for multipartite 
entanglement: the hopping term, which spreads the wavefunction [a state with on-site localization obviously shows no 
global entanglement], and the interaction term, which further enhances the correlations between the particles. High 
levels of entanglement require more than simply derealization. As we discussed in Sec. II. E, it is possible to have 
very delocalized states with low entanglement. 

Also in contrast with NPC, when A < 1, global entanglement decreases monotonically with disorder. This is a result 
of the breaking of symmetries. It is possible, although this requires further investigation, that besides S 2 symmetry 
(and consequent rotational symmetry), reflection also has positive effects on quantum correlations. To illustrate this 
idea, compare the values of NPC and Q in the table below. We have a system with L = 12, A = 0.5, and two equal 
defects at the borders, u>i — lol = 0.05 J, which guarantees that the S 2 symmetry is broken, but not parity. Suppose 
that we may include an additional defect on site n with energy 0.5J (if n = L, then uj^ becomes 0.55J). This is a 
simple way of breaking parity. 





NPC 


Q 


No additional defect 


217 


0.989 


Additional defect at n = L 


176 


0.953 


Additional defect at n = (L/2) + 1 


256 


0.987 



The chain with equal or unequal border defects is integrable [67|, whereas a defect in the middle leads to chaos [251 ] . 
Independently of the position of the additional defect, it simply decreases Q. NPC, on the other hand, decreases 
only in the integrable regime, but increases in the chaotic region. It has been argued in Ref. [28[ that Q seems closer 
related to NPC than with the integrable-chaos transition. Here, we go one step further and claim that, even though 
derealization is necessary for the existence of quantum correlations, the two quantities do not always go hand in 
hand. The onset of chaos has a positive effect on spatial derealization, whereas for global entanglement the breaking 
of symmetries associated with chaos are more detrimental than possible gains associated with further derealization. 

Large interactions (A > 1) are unfavorable to both NPC and Q, since the spectrum becomes gapped. However, 
for anisotropies not too large, A > 1, as d/J increases and the energy bands overlap, the Q curves quickly surpass 
the ones for smaller interactions. This reinforces once again the fundamental role of interactions in the creation of 
quantum correlations. 

2. FP- and IP-basis 

The results for NPC are entirely dependent on the basis chosen for the analysis. While the site basis is physically 
motivated by studies of spatial localization, further insight may be gained by considering alternative bases. To 
separate the regular motion dictated by integrable Hamiltonians from chaoticity, we show in Fig. [3] the behavior of 
derealization vs. disorder for NPC computed in the FP- and IP-basis. 
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FIG. 3: Average NPC vs. uncorrelated disorder. Left: NPC in the free-particle basis; from bottom to top: A = 0, 0.5, 1, 1.5, 2, 5 
Right: NPC in the basis of the integrable interacting clean chain. Average over 924 states and 20 realizations; L — 12, M = L/2; 
a = 0. 



The clean one-dimensional system considered here, with or without interactions, is integrable, chaos emerges only 
when disorder is added to it. In the FP-basis, a possible correlation between chaos and the complexity of the 
wavefunctions is reflected in the rate of derealization, which increases abruptly in the chaotic region of < d/ J < 0.5, 
< A < 1, and is less dramatic for A = and A > 2. But overall, NPC simply increases with A. In contrast, 
in the IP-basis the behavior of NPC with A for < d/J < 0.5 is non-monotonic. Here, as in the site-basis case, 
derealization is maximum where rj is minimum, the largest value appearing again for A = 0.5. It is only at large 
values of disorder that the behavior of NPC with A becomes trivial. Therefore, the effects of the interplay between 
interaction and disorder and the consequent onset of chaos are well singled out by the non-trivial behavior of NPC in 
both site- and IP-basis. 



B. Long-range correlated random disorder 



We now focus on a disordered chain with d/J — 0.25 - the chaotic region of interacting systems - and proceed with 
a comparison between uncorrelated and correlated random disorder. 
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FIG. 4: Interplay between interaction and long-range correlated random disorder. Left panel: Average derealization in the 
site-basis vs. a. Right panel: Average global entanglement vs. a. All panels: d/J — 0.25, L = 12, M = L/2, average over 924 
states and 20 realizations. 



The left panel of Fig. 5] shows the average NPC vs. a. The potential for spreading the wave functions associated 
with a > has a non-monotonic behavior with the anisotropy. This may be explained as follows. Correlated disorder 
indicates more order, in a sense, it brings the system closer to a clean chain; but how the system approaches its clean 
limit is highly dependent on the interaction amplitude and the value of NPC we start with when a = 0. Long-range 
correlated disorder can always increase derealization when A = 0, since a system with free particles is maximally 
delocalized in the absence of disorder (main panel in Fig. [T]). The growth in NPC is obviously not sufficient to reach 
the values of a clean system, and not even those obtained for chaotic systems with A < 1, but, as shown in Fig. [4] 
it approaches the level of derealization of a chaotic chain with A ~ 1. For interacting systems, on the other hand, 
correlated disorder may decrease NPC. The maximum value of NPC when a = depends on the amplitude of the 
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uncorrelated disorder, d/J. When d/J — 0.25 and a = 0, it is seen from the main panel in Fig. [T] that the NPC 
curves with A < 1 have already passed their maximum, while those with A > 1, such as A = 2, are right at their 
peak. Therefore, for A < 1, by increasing a we can further delocalize the system, whereas for A > 1, correlated 
disorder has the opposite effect of decreasing NPC. For large values of A, such as A = 5, since NPC simply decays 
with uncorrelated disorder (main panel of Fig. [1]), correlated disorder will always increase derealization. 

Insight on the non-trivial behavior of NPC with a may be gained by studying, in Fig. [H how correlated disorder 
affects the histograms for the diagonal elements and their consequences on the relation between NPC and the 
eigenvalues. We select four values of the interaction amplitude, A = 0, 0.5, 1.3, and 2.0. 




FIG. 5: Uncorrelated vs. long-range correlated random disorder. First and second rows: normalized histograms for the diagonal 
matrix elements Hoi with uncorrelated and correlated disorder, respectively. Third row: NPC vs. eigenvalues; black curve: 
a = 0; red curve: a = 10. Columns from left to right: A = 0, 0.5, 1.3, 2. All panels: L = 12, M = L/2, 40 realizations. 

For uncorrelated disorder, as A increases, the distribution of simply broadens (Fig. [5J ai, bi, Ci, di), however 
NPC does not decrease monotonically. As seen in the main panel Fig. [IJ for d/J — 0.25, A = and 1.3 lead 
approximately to the same value of NPC, derealization is maximum at A ~ 0.5 and is very low for A = 2. This 
non-trivial behavior is associated with the transition to chaos from A = to A = 0.5: The featureless curve of NPC 
vs. E (as), which shows intermediate values of derealization for all energies, is substituted by a curve that peaks in 
the middle of the spectrum (b^). As expected from chaotic TBRE's, for A = 0.5, NPC approaches the GOE value 
of N/3 in the middle of the spectrum and has smaller values only at the edges. This results in an overall increase of 
(NPC). It is only from A = 0.5 to A = 1.3, 2.0, that the broadening of the distributions of diagonal energies leads to 
the decrease of NPC (03^3). 

The inclusion of random disorder with long-range correlation (a — 10) has little effect on the derealization of 
a half- filled system with interacting particles, but it significantly affects the model with non-interacting particles. 
When A = 0, correlated disorder increases the number of resonances in the middle of the histogram (compare ai and 
a2), leading to an interesting shaped distribution. This is reflected in the larger values of NPC in panel a3. When 
A = 0.5, 1.3, the shape of the distribution changes slightly (cf. bi, b2 and Ci, C2), having little effect on the curves 
for NPC vs. E (b3, C3). However, when A = 2 the energy band overlapping obtained by adding uncorrelated disorder 
(see a;, bi, and c; of Fig. [1] and di here) is partially removed by adding order via a (d2). This explains the decay 
of (NPC) verified in Fig. 0] The effects of correlated disorder are, however, limited and all curves for NPC and Q 
saturate for a > 4 (see Fig. 2]). 

There are then two possibilities to increase derealization in a chain of free particles with uncorrelated disorder: 
by including weak interaction, so that chaos may set in, or by pushing the system toward the clean limit by adding 
correlated on-site energies. The two alternatives constitute, however, opposite procedures, so when put together their 
effects do not add up. 

In general, the behavior of multipartite entanglement and derealization with correlated disorder are comparable. 
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Similarly to NPC, the effects of a on Q also depend on the value of Q we start with when a = 0, an information 
extracted from Fig. [5] In the presence of correlated on-site energies, the right panel of Fig. 0] shows that (i) global 
entanglement increases with a when A = 0, (ii) Q is not much affected by correlated disorder when A ~ 1, (iii) Q 
decreases with a when A > 1. Notice, however, the prominent role of interaction in establishing quantum correlations. 
In particular, compare the behavior of the curves for A = and 1.3. When a = 0, NPCa=o ~ NPCa=i.3, while 
Qa=i.3 ~ 1-06 Qa=q- When a = 8, NPCa=o increases significantly and becomes approximately 20% larger than 
NPCa=i.3, while the growth of Qa=o is more limited and it remains smaller than the global entanglement for A = 1.3, 
<2a=i.3 ~ 1-03 <3a=o- 



IV. DILUTE LIMIT 



The dilute limit implies M <C L. Here, we focus on the smallest value of M where the interaction plays a role, 
the two particle case, M = 2, and study the competition between interaction and both uncorrelated and correlated 
disorder. 




FIG. 6: Top panels: interplay between interaction and uncorrelated disorder, a = 0; bottom panels: interplay between 
interaction and long-range correlated disorder, d/J — 0.25. Top left (right) panel: Derealization (Global entanglement) vs. 
anisotropy for different values of d/J. Bottom left (right) panel: Derealization (Global entanglement) vs. a for different values 
of A. Left bottom panel - curves from top to bottom at a = 8: A = 0, 0.25, 0.5, 0.75, 1, 1.35, 2. All panels: L = 24, M = 2, 
average over all states and 20 realizations. 



From the top panels of Fig. [51 one sees that for a fixed value of the interaction amplitude, NPC and Q decay 
monotonically with disorder. This justifies our choice here to consider A in the abscissa, instead of d/J, as in Figs.[TJ 
[2J For d/J^0, such as d/J = 0.25 and 0.5, a value of A ^ may still exist where NPC and Q become maximum, but 
the effects of the interaction are now much less prominent than in the half-filled chain. In the dilute limit, the H\y 
term dominates over the Ising interaction and becomes the main determining factor for the level of derealization and 
global entanglement. As a result, the behaviors of NPC and Q are more comparable than at half-filling (cf. Fig. [6] 
and Figs. [HE]). 

The bottom panels show how NPC and Q behave with correlated disorder for various values of A, keeping fixed 
d/J — 0.25. Contrary to the half-filled case, for all cases considered here, NPC and Q increase with a. This is a 
reflection of the trivial behavior of NPC and Q with uncorrelated disorder. Moreover, the effects of a are now more 
pronounced than in the half-filled chain, especially for the chain with A = 0. NPCa=o now surpasses all curves shown 
and even Qa=o manages to outperform Qa=i.35- In comparison to the M = L/2 case, correlated disorder can now 
bring NPC and Q closer to the values of a clean chain. These results simply reflect the minor role of the interaction 
in a dilute system. The increase of NPCa=o from a = to a > 4 is more than 50% for L = 24, M = 2. As expected, 
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this behavior is yet amplified in larger systems, as shown in the left bottom panel of Fig. [7] 
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FIG. 7: Effects of the chain size on derealization in the presence of long-range correlated disorder. Top panels: M = 1, bottom 
panels: M — 2 and A = 1. Left panels: Derealization vs. a for different values of L; right panels: Derealization vs. chain 
size for different values of a, curves from bottom to top: a = 0, 1, 2, 3, 4, 5. Average over all states and 40 realizations. 



The right panels of Fig. [7] compare derealization with system size in the case of M = 1 and M — 2 for different 
values of a. When M = 1 the dimension of the Hamiltonian is N = L. For a = 0, as seen in the top right panel, 
apart from small system sizes where border effects are important, there is little difference between the NPC values 
for different L's. This result is in agreement with the fact that Anderson localization occurs in one-dimensional 
systems with a single particle, that is, uncorrelated random on-site energies lead to a finite localization length of the 
wavefunctions. On the other hand, for large values of a, the increase rate of NPC becomes linear. Large correlations 
between on-site energies bring the behavior of NPC closer to that of a clean system, where NPC oc L. This is consistent 
with Ref. [8j], where it was shown that, in the thermodynamic limit, one-particle states become delocalized for a > 2. 

When M = 2, the dimension of the Hamiltonian is quadratic on the system size, N = L(L — l)/2. In this case, as 
shown in the right bottom panel, NPC grows quadratically with L even when a = 0. This suggests that, on average 
and for the finite systems studied here, there is no sign of Anderson localization for states with two interacting 
particles, even when the disorder is uncorrelated. There are, however, few isolated states with small NPC (data not 
shown), which accounts for situations where the particles are sufficiently apart, so that the interaction is negligible 
and the particles are individually localized. Further studies are now required to conclude if, on average, localized 
states will be the dominating scenario for larger chains, or if interaction will in fact prevent localization even in the 
thermodynamic limit. 



V. DISCUSSION AND CONCLUSION 



We studied how derealization and global entanglement in a quantum many-body one-dimensional system are 
affected by the interplay between on-site (uncorrelated and long-range correlated) random disorder and nearest- 
neighbor interaction. The analysis applies equally to a chain with spins-1/2, spinless fermions, or hard-core bosons. 
Half-filled and dilute chains were considered. 

In a half-filled chain, the largest value of spatial derealization (NPC) appears in a clean chain (d/J = 0) with 
non-interacting particles (A = 0). Here, the only constraints in the spreading of the wavefunctions are the symmetries 
of the system and the boundary conditions. The addition of uncorrelated disorder to a non-interacting system or the 
inclusion of interaction to a clean chain, suppresses resonances and decreases NPC. On the other hand, when both A 
and d/J are present, the behavior of NPC becomes non-trivial. A clean chain with interacting particles constitutes 
an integrable model with eigenvectors delocalized in the site-basis. As uncorrelated disorder is added, symmetries are 
broken. In the case of weak interaction, < A < 1, this leads to the onset of chaos and the further increase of NPC 
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in the site-basis. 

The behavior of global entanglement (Q) is influenced by derealization, but does not follow it exactly. In particular, 
for d/J — ► 0, the largest values appear in the presence of weak interaction. Also different, in weakly interacting systems, 
is its steady decay with uncorrelated disorder and the consequent breaking of symmetries. Besides derealization, key 
contributing factors for global entanglement are then interactions and symmetries. 

We have also briefly commented on the results for NPC in bases other than the site-basis. The latter is considered 
in studies of spatial localization, which was our main goal in this work. However, it is also useful to separate regular 
motion from chaoticity by analyzing derealization in a basis where the model is integrable. In the basis consisting of 
eigenstates of a clean system with interaction (IP-basis), a non-monotonic behavior of NPC with A was also verified 
for disordered systems with d/J < 1. This non-trivial behavior of NPC with interaction in both the IP- and the 
site-basis appears to be a main indication of the complexity enhancement of wavefunctions in the chaotic region. 

The presence of correlated disorder (a ^ 0) has different consequences depending on the value of A. By correlating 
on-site energies, the level of order increases. A disordered system with non-interacting particles then gets closer to 
the limit of a clean chain and NPC increases. In the presence of interaction, the effect is opposite, since here it is 
disorder that enhances derealization by breaking symmetries and overlapping energy bands. The behaviors of NPC 
and Q with a reflect the non-trivial behavior of these quantities with d/J when a = 0. 

The study of a dilute system allows for the consideration of larger chains. In the presence of uncorrelated disorder, 
contrary to the one-particle case, on average, our results show no sign of Anderson localization in the finite systems 
considered. Long-range correlation delocalize both one and two-particle states. 

In a future work, we intend to search for conclusive evidences for the absence (or not) of localization of two- 
interacting-particle states in systems with uncorrelated disorder. We also aim to study in details the effects of 
symmetries on non-local correlations. Another question worth investigation is how the integrable and chaotic regimes 
of the considered one-dimensional systems relate to different transport behaviors, namely diffusive or ballistic 68], 
and to the problem of thermalization [HJ, [69| . Also interesting would be the comparison of results with the Hubbard 
model and the inclusion of long-range interactions to the systems considered. 

We hope that our theoretical results will motivate experimental verifications. Useful tools to simulate many-body 
effects in simplified models of condensed matter physics as the ones described here are optical lattices. They allow 
for control of the parameters of the system, such as the strength of interaction and the level of disorder [7(| • 
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